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I.  INTRODUCTION 


In  the  present  report  we  study  the  velocity  end  temperature  boundary-layer 
development  within  a  centered  expansion  wave  moving  into  a  stationary  fluid.  A 
centered  expansion  wave  can  be  generated  either  in  a  shock  tube  or  e  tube  wind 
tunnel  such  as  designed  by  Ludwig. *  It  is  assumed  that  the  fluid  has  initially 
a  uniform  temperature  TQ  but  is  separated  by  a  diaphragm  in  two  regions,  with 
different  pressures.  If  the  diaphragm  is  suddenly  removed,  an  expansion  wave 
propagates  into  the  stationary  high  pressure  region.  In  the  expansion  region, 
the  temperature  of  the  fluid  decreases  due  to  the  drop  in  the  pressure. 
Therefore,  if  the  wall  is  kept  at  the  initial  temperature  TQ,  there  will  be  a 
net  heat  transfer  from  the  wall  to  the  fluid.  This  heat  transfer  and  the 
boundary-layer  development  in  the  expansion  region  is  the  subject  of  the  present 
analysis. 

2 

As  early  as  1859,  Riemann  analysed  the  propagation  of  a  plane  disturbance 
of  finite  amplitude  in  an  unsteady  one-dimensional  isentropic  flow. 

Subsequently, many  other  authors  studied  the  shock  tube  problem.  Huber  et  al , 5 
for  example,  provided  analytical  expressions  for  the  velocity,  pressure,  and 
temperature  of  the  inviscld  flow.  However,  in  order  to  calculate  the  heat 
transfer  characteristics,  the  viscous  boundary- layer  development  on  the  inner 
wall  must  be  computed.  Mirels^  analyzed  the  boundary-layer  flow  behind  the 
expansion  wave  by  assuming  that  the  expansion  wave  is  a  line  of  discontinuity 
separating  the  undisturbed  high  pressure  region  from  the  rarefied  gas  region. 
Cohen5  used  a  coordinate  expansion  method  to  solve  the  boundary-layer  flow 
within  the  centered  expansion  wave  with  a  finite  width.  Cohen's  treatment, 
which  is  an  Improvement  over  Mirels'  solution,  is  based  on  the  assumption  that 
the  wall  is  isothermal  or,  equivalently,  that  the  thermal  conductivity  of  the 


^H. Ludwig,  "Tube  Wind  Tunnel:  A  Special  Type  of  Blowdown  Tunnel, "  AGARD 
Report  143 ,  July  1957. 

2B.  Riemann,  "Ueber  die  Fortpflanzung  ebener  Luftwellen  von  endlicher 
Schwingungeweite ,"  Abh.  dev  Koniglichen  Geseellechaft  dev  Wieeenechaften , 
Vol.  8,  pp.  43-65,  1858-1859 . 

3P.W.  Huber,  C.E.  Fitton,  and  F.  Delpino,  "Experimental  Investigation  of 
Moving  Pressure  Disturbance  and  Shock  Waves  and  Correlation  wiih  One- 
Dimensional  Unsteady-Flow  Theory,"  NACA  TH  1903,  1949. 

4H .  Mirels,  "Boundary  Layer  Behind  Shock  or  Thin  Expansion  Wave  Moving  into 
Stationary  Fluid,"  HACA  TN  3712,  1956. 

**N.B.  Cohen,  "A  Power-Series  Solution  for  die  Unsteady  Laminar  Boundary-Layer 
Flow  in  an  Expansion  Wave  of  Finite  Width  Moving  Through  a  Gas  Initially  at 
Rest,  "  SAC  A  TH  3943,  1957. 
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vail  la  infinite.  Hall*  later  generalised  Cohen's  solution  by  allowing  the  wall 
to  have  a  finite  theraal  conductivity,  and  gave  solutions  of  heat  transfer  for 
SOM  representative  values  of  theraal  conductivity.  Hall's  solution  in  the  case 
of  adiabatic  wall  predicts  a  teaperature  in  the  boundary  layer  that  is  equal  to 
the  teaperature  in  the  freeatreaa. 

The  present  study  provides  solutions  for  the  boundary-layer  flow 
corresponding  to  two  wall  theraal  conditions,  namely,  an  adiabatic  condition  and 
am  isothermal  condition  at  undisturbed  fluid  temperature.  The  solutions  are 
obtained  by  applying  Howarth's  transformation  and  a  similarity  transformation 
via  one-paraMter  groups  which  reduces  the  number  of  Independent  variables  from 
three  to  two  and  then  a  series  expansion.  The  first  three  terms  of  the  series 
are  used  to  obtain  a  solution  of  second-order  accuracy. 


II.  PROBIEM  FORMULATION 

The  boundary-layer  flow  in  the  region  of  the  expansion  wave  as  shown  in 
Figure  1  la  considered  to  be  two-dimensional,  unsteady,  compressible,  and 
laminar.  The  fluid  is  ideal  gas.  External  to  the  boundary  layer,  the  flow  is 
assumed  to  be  unsteady,  one-dimensional  and  invlscid.  The  freestream  flow  is 
assumed  to  be  known  and  it  is  used  as  the  matching  condition  for  the  boundary- 
layer  solution. 


6 

J .G.  Hall ,  . Laminar  Boundary  Layer  Developed  Within  Unsteady  Expansion  and 

Compression  Waves, "  AIAA  Journal,  Vol.  10,  p.  499,  April  1972. 

?L.  Howarth,  "Concerning  t he  Effect  of  Compressibility  of  Laminar  Boundary 
Layers  and  Their  Separation,  "  Proceedings  of  the  Royal  Society,  Series  A, 
Vol .  194,  p.  16,  1948. 

O 

A.J.A.  Morgan,  "The  Reduction  by  One  of  the  Number  of  Independent  Variables 
in  Some  System  of  Partial  Differential  Equations,"  Q.J.M.,  p.  250 ,  1952. 

9 

A.G.  Hansen,  Similarity  Analyses  of  Boundary  Value  Problems  in  Engineering, 
Prentice  Hal l,  Englewood  tfliffs,  NJ,  1964. 
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Figure  1.  Shock  Tube  Flow  at  t  -  t* 


We  assume  that  the  undisturbed  region  is  kept  at  a  temperature  TQ.  Then 
the  wavefront  of  the  expansion  wave  propagates  at  the  sound  speed  of 
Ur  -  (yRTq)^^.  The  velocity,  temperature,  and  pressure  in  the  expansion  region 
outside  the  boundary  layer  are  given  by  Huber  et  al  ^  in  terms  of  a  similarity 


variable  5  as  follows: 


Velocity 


ui<t’lt)  "  r+T  5 


F(5) 


Temperature 


o^t.x) 


—  ( 1  -  III  r  y 


G(5) 


Pressure 


PjCt,*) 


Y_,  2y/ (y-1) 

-■<«> 


(I) 


(2) 


(3) 


where  the  similarity  variable  £,  derived  from  the  similarity  transformation  via 
one'pararaeter  groups  is 


C  -  1  +  (x/t) 


x  -  x/L 


(4) 


and  L  la  a  reference  length.  As  shown  In  Figure  1,  £-0  denotes  the  expansion 
wavefront  and  €-1  Is  located  at  the  origin  of  the  wave.  F(£).  G(£),  and  H(£) 
are,  respectively,  called  the  velocity  function,  temperature  function,  and 


pressure  function  of  the  Invlscid  flow.  The  velocity  u^  is  the  inviscid 
velocity  normalized  by  Ur,  the  temperature  0. by  the  reference  temperature 

T  _  T]2  /  d  .  II*  n  1  -  l  1 1 »  m «-  .  * .  . 


rr  -  Uy/R,  and  the  pressure  p^  by  PQU^. 


is  the  density  of  the  gas  at  rest 


ahead  of  the  expansion  wave. 


In  the  boundary  layer  we  assume  that  the  thermal  conductivity  k  and 
specific  heat  C  are  constants  and  that  the  viscosity  y  is  proportional  to  the 
temperature.  Tne  viscosity  y  can  be  made  dimensionless  as 


y  -  y/yr  -  T/Tr 


(5) 


where  the  reference  viscosity  yf  is  evaluated  at  the  reference  temperature  Tr. 
We  also  define  the  following  dimensionless  variables  and  parameters: 

1  /2 

y'  -  Re  '  Y/L,  u  -  U/Ur  , 

1/2  2 

v  “  Re  V/Ur,  p  -  P/(pQUr  ),  0  -  T/Tr,  p  -  p/pQ  ,  (6) 

Re  -  UrL  P<j/yr,  Pr  -  Cpyr/k,  Tr  -  U^/R  . 

The  governing  equations  for  the  two-dimensional,  unsteady,  compressible 
boundary-layer  flow  are,  in  terms  of  these  variables. 

Continuity 

+  h  (pu)  +  ip* (pv) " 0  <7> 

Momentum 


_  i£  +  _L_ 

dx  3y  * 


(u 


(8) 


10 


State 


P  ■  P0 


(10) 


The  pressure  p  In  Eq.  (8)  is  the  same  as  the  p^  in  Eq.  (3). 
Howarth's  transformation7 


ter  introducing 


y  -  /  pdy ' 
o 


(ID 


a  streamfunction  i|>,  which  satisfies  the  continuity  equation,  can  be  defined  as 


u  3y,v  pl3x+3y3x  +  3t^ 


(12) 


Substituting  Eqs.  (11)  and  (12)  into  Eqs.  (8)  and  (9),  we  have 
atay  3y  3x3y  3x  3y2  p  dx  p  3y3 


and 


+  ,  i^i  ,  e* 

3t  3y  3x  dx  3y  "  Pr  2  y  '•p' 

®y 


r3p  +  3+  3p  p^  32»  2 

lat  ay  3x  6  \  2}  J 

ay 


(13) 


(14) 


In  addition,  we  have  the  following  boundary  conditions  for  and  0: 


At  the  wall,  y  ■  0 


u  ■  ■  0  (no  slip  condition) 

v  «  — -  4^  “  0  (impermeable  condition) 
p  3x 

0-0  -  —  (isothermal  wall  at  T  ) 

w  y  o 

or 

-|y-  «  0  (adiabatic  wall) 

At  the  outer  edge  of  the  boundary  layer,  y  >  « 

u  ■!£- 
0  «  01(t,x) 


(15) 


(16) 


(17) 


(18) 

(19) 


III.  METHOD  OF  SOLUTION 

We  obtain  a  solution  of  Eqs.  (13)  and  (14)  through  a  similarity 
transformation  via  one— parameter  groups  and  a  subsequent  series  expansion  of  the 
functions  <J>  and  0.  The  similarity  transformation  reduces  the  three  independent 
variables  (t,x,y)  to  two  similarity  varigbles  (£,ri)  which  are 

£  ■  1  +  (x/t)  and  n  *  y/t*^  (20) 


A  detailed  derivation  of  these  two  similarity  variables  is  given  by  Chang. 

The  similarity  transformation  also  transforms  the  dependent  variables,  ij>  and  0 
in  Eqs.  (13)  and  (14)  into 


10L.M.  Chcmg,  "Unsteady  Compressible  Boundary  Layer  Flows  Within  a  Centered 
Expansion  Wave  and  Behind  a  Strong  Plane  Blast  Wave,"  Ph*D»  Dissertation, 
The  University  of  Iowa ,  December  1972 . 
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(21) 


► , 


*  % 


♦<t,x,y>  ■  t1/2f({,n) 


I 

>  » 

L* 

\ 


i 

p* 

r. 

I\ 

I 


I 


I 


r- 


e(t,x,y)  »  g(5.n) 


(22) 


Substitution  of  Eqs.  (3)  and  (20-22)  into  Eqs.  (13)  and  (14)  results  in 


(l-£)fr  -  ^  nf  +  f  fr  -  f_  f  ■  -  g  +  Hf 

2  nn  n  Sn  5  nn  H  nnn 


(23) 


and 


u-c)*5  -j  +  fng5 


f5*n  "  Pr  Hgtm 


(24) 


+  ^  «  [(l-C)H'  +  f  H*  +  —  f2  J 
Y  H  L'  n  g  nnJ 


where  H'  is  the  derivative  of  H(£)  with  respect  to  C*  A  repeated  application  of 
the  transformation  technique  shows  that  no  new  similarity  variables  that  combine 
the  variables  E,  and  n  exist.  However,  we  can  define  a  locally  similar  variable 
as 


z  “  n/£ 


1/2 


(25) 


Using  the  new  variable  z,  a  fast-converging  series  solution  to  the  problem  may 
be  constructed  as  follows 


f  -  51/2F(U  *  l  £nf  (z) 


nm0 


■-£-£-««>  j,,  «-«;<■> 


a Q 

e  -  g(5,n)  -  7  +  [g(0  -  -  -  l  5ngn(z) 

n»0 


(26) 


(27) 


(28) 


HL'VV.x'fv1' 


Substituting  the  above  expressions  and  their  derivatives  into  Eqs.  (23)  and 
(24),  we  have,  after  grouping  the  terms  with  the  same  power  of  £,  the  following 
system  of  ordinary  differential  equations. 
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rr . 


Zero-order  equations  (n»0) 


fjj*  +  (l/2)zf"o  -  f*  -  -  1 


(1/Pr)g"  +  ( 1/2) zg'  -  g  -  -  1 


First-order  equations  (n»l) 


f”'+—  m 

*12  Ztl  "l  Y+l  *0 


yTT  fofo  -  t'o 


+  —2—  ( f  *  >  ^  +  2  ni  g  -  XzJL 

Y+l  '“o'  Y+l  go  Y+l 


1 _  *«  + 1  ZK<  .  2k  -  2-  -ii-A-  *-  +  ri  nil  K* 

Pr  gl  2  gl  gl  Pr  2(y+1)  go  U  y+1J  go 


-  (^-)  sD  -  (^r)  8of0  +  (4r)  fo 


-  (^r)»;  +  (4r)f>o + 


Second-order  equations  (n-2) 


f;*  +  ( 1/2) zf 2  -  3f 2  -  F*(z) 


(1/Pr)g”  +  (l/2)zg^  -  3g2  -  G*(z)  (34) 

The  expressions  F*(z)  and  G*(z)  are  given  in  the  Appendix.  The  boundary 
conditions ,Eqs.  (15)  -  (19), become: 


at  the  wall,  z  *  0,  for  n  >  0 


fn<°)  -  0,  f;(0)  -  0 


gnCO)  -  0 


for  Isothermal  wall 


(35) 


or 


g;(0)  ■  0  for  adiabatic  wall 

at  the  outer  edge  of  the  boundary  layer,  z ■+**> 

f»  («»)  >  i  f'(«»)  -  0  for  n>l 

o  n 

g0(«)  ■  1.  gn(“)  “0  for  n>l  (36) 

g'(«o)  -  0  for  n>0 

n 


IV.  NUMERICAL  RESULTS  AND  DISCUSSION 

The  preceding  three  sets  of  ordinary  differential  equations  were  solved 
numerically.*®  With  Pr  «  0.72  and  y  *  1.4  the  numerical  results  of  the  velocity 
functions  f^,  fj,  and  f£,  and  the  temperature  functions  gQ,  gj,  and  g2  for  both 
isothermal  and  adiabatic  walls  are  tabulated  in  Tables  1  and  2  and  plotted  in 
Figures  2  and  3.  Since  there  is  no  temperature  term  involved  in  the  zero-order 
momentum  equation,  Eq.  (29),  the  function  f^  is  independent  of  the  temperature 
condition  at  the  wall.  Therefore,  as  shown  in  Figure  2,  the  velocity  function 
f^  is  the  same  for  the  two  wall  conditions.  However,  the  function  fj  and  f £  are 
certainly  affected  by  the  wall  temperature.  For  the  case  of  isothermal  wall  the 
present  results  are  Identical  to  Cohen's  solutions.*® 

Once  the  zero-order,  first-order,  and  second-order  functions  are  computed, 
the  solutions  to  the  problem  in  a  three-term  series  can  be  readily  constructed. 


TABLE  1.  VELOCITY  FUNCTIONS 


z 

Isothermal  wall 

Adiabatic  wall 

fo 

/; 

t'i 

f'o 

/; 

/] 

0 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.2 

0.20643 

0.13364 

0.06425 

0.20643 

0.11252 

0.05276 

0.4 

0.37735 

0.22220 

0.09358 

0.37735 

0.18962 

0.08014 

0.6 

0.51717 

0.27376 

0.09495 

0.51717 

0.23642 

0.08529 

0.8 

0.63011 

0.29602 

0.07631 

0.63011 

0.25834 

0.07322 

1.0 

0.72016 

0.29611 

0.04555 

0.72016 

0.26064 

0.04973 

1.2 

0.79100 

0.28046 

0.00973 

0.79100 

0.24907 

0.02052 

1.4 

0.84596 

0.25460 

-0.02540 

0.84596 

0.22773 

-0.00945 

1.6 

0.88801 

0.22310 

-0.05576 

0.83800 

0.20081 

-0.03638 

1.8 

0.91971 

0.18957 

-0.07788 

0.91971 

0.17156 

-0.05778 

2.0 

0.94325 

0.15666 

-0.09373 

0.94325 

0.14245 

-0.07243 

2.4 

0.97288 

0.09917 

-0.10045 

0.97288 

0.09087 

-0.08185 

2.8 

0  98783 

0.05719 

-0.08558 

0.98783 

0.0S27I 

-0.07162 

3.0 

0.99205 

0.04202 

-0.07417 

0.99205 

0.03882 

-0.06261 

3  4 

0.99680 

0.02128 

-0.05023 

0.99680 

0.01975 

-0.04287 

3.8 

0.99848 

0.00989 

-0.03016 

0.99884 

0.00925 

-0.02585 

40 

0.99935 

0.00653 

-0.02246 

0.99935 

0.00615 

-0.01924 

5  0 

1.00010 

0.00047 

-0.00389 

1.00010 

0.00062 

-0.00317 

6.0 

1  00020 

0.00018 

-0.00696 

7.0 

1.00020 

-0.00064 

-0.001 10 

8  0 

1  00020 

-0  00099 

-0.00125 

TABLE  2.  TEMPERATURE  FUNCTIONS 


Isothermal  wall 

Adiabatic  wall 

< 

fa 

fi  fa 

fa 

ft 

« 1 

0 

0.00000 

0.00000 

0.00000 

1.00000 

-0.18005 

0.2 

0.17756 

003701 

-042294 

140000 

-0.17984 

0.4 

0.32907 

047981 

-041369 

1.00000 

-0.16077 

0.6 

0.45723 

0.11981 

0.00384 

140000 

-0.13707 

0.6 

0.56460 

0.15307 

0.01733 

140000 

-0.11271 

1.0 

0.65394 

0.17434 

0.02195 

1.00000 

-0.09000 

1.2 

0.72739 

0.10629 

0.01737 

1.00000 

-0.07014 

1.4 

0.70705 

0.10800 

0.00563 

1.00000 

-0.05352 

1.6 

0.13554 

0.18346 

-0.01035 

1.00000 

-0.04009 

1.0 

0.87411 

0.17213 

-042768 

1.00000 

-0.02954 

2.0 

0.90460 

0.15673 

-0.04395 

1.00000 

-0.02144 

2.4 

0.94689 

0.12032 

-0.06726 

1.00000 

-0.01004 

2.1 

0.97167 

0.00461 

-0.07489 

1.00000 

-0.00522 

34 

0.97969 

0.06091 

-0.07341 

1.00000 

-0.00356 

3.4 

0.90994 

0.04334 

-0.06322 

1.00000 

-0.00159 

3.1 

0.99532 

0.02553 

-0.04068 

1.00000 

-0.00068 

4.0 

0.99691 

0.01916 

-0.04068 

1.00000 

-0.00436 

S.0 

0.99996 

0.00383 

-0.01241 

1. 00000 

6.0  1.00050  0.00011 

0.00033 
0.00007 


-0.12156 

-0.10631 

-0.07399 

-0.03025 

-0.00605 


0.03229 

0.04042 

0.04205 

0.04134 

0.03744 

0.02695 

0.01717 

0.01323 

0.00742 

0.00309 

0.00275 

0.00041 


-0.00270 

-0.0010! 

-0.00100 
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Figure  2.  Profiles  of  Velocity  Functions 


*  °  #  *  %.9i!&9*  #  ‘  M 


Figure  3.  Profiles  of  Temperature  Functions 


A.  Velocity  Profiles 

From  Eq.  (27),  the  velocity  in  the  boundary  layer  is  simply 


u  -  ¥{f'0  +  Cf*  + 


u/ul  "  fo  +  5fi  +  ?2f2  (38) 

for  F  ■  u^  given  in  Eq.  (1).  Figures  4  and  5  show  u/u^  vs  z,  respectively,  for 
Isothermal  and  adiabatic  walls.  One  observes  that  at  a  given  z,  u/u^  increases 


with  (.  This  Is  because,  as  shown  in  Eq.  (1),  the  invlscld  velocity  u^  relative 

to  the  wavefront  appears  to  accelerate  in  the  positive  £  direction.  In  other 

words,  the  fluid  is  expanding  from  the  high  pressure  to  the  lower  pressure 

region  and  the  further  the  downstream  distance  is  from  the  wavefront  the  higher 

the  fluid  velocity  becomes.  As  a  result,  even  though  the  viscosity  diffusion  in 

the  boundary  layer  tends  to  increase  the  boundary  growth  the  acceleration  may 

suppress  some  growth  of  the  boundary  layer.  If  the  boundary-layer  thickness  is 

taken  to  be  approximately  z.  ■  4  from  Figures  4  and  5,  Eqs.  (25)  and  (20)  give 

6 


*5  '  •  4 


(39) 


Thus, at  a  given  instance,  the  boundary-layer  thickness  behind  the  wavefront 
may  grow  like  4/x+t.  The  larger  the  time,  the  slower  the  growth  of  the 
boundary-layer  thickness  will  be  with  respect  to  x.  In  other  words,  the 
boundary  layer  grows  rapidly  immediately  behind  the  expansion  wavefront  and 
flattens  out  further  downstream  from  the  wavefront. 


Uflli 


Figure  4.  Velocity  Profiles  for  an  Isothermal  Wall  at 
£  -  0,  0.1,  0.2,  and  0.3 


Figure  5.  Velocity  Profiles  for  an  Adiabatic  Wall  at 
£  ■  0,  0.1,  0.2,  and  0.3 


A  careful  comparison  between  the  profiles  shown  in  Figures  4  and  5  reveals 
that  at  a  given  s  the  velocity  is  slightly  higher  for  the  isothermal  wall  due  to 
heat  addition  from  the  wall  to  the  fluid. 

B.  Temperature  Profiles 

The  temperature  in  the  boundary  layer  is 

e  -  7  +  (g  -  i)  (g0  +  5gl  +S2g2) 

or 

f--  [7+  (G  (gD  +  5g A  +  C2g2)]/G  (40) 


The  solutions  are  plotted  in  Figure  6.  The  temperature  in  the  boundary  layer  is 
seen  to  be  higher  than  that  in  the  freestream  and  the  value  of  6/6^  increases 
with  5  as  a  result  of  an  Increasing  viscous  heating  in  the  boundary  layer. 

Figure  6  shows  that  at  £  *»  0.3  the  wall  temperature  is  higher  than  the  local 
lnvi8cid  temperature  by  approximately  10  percent. 
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Figure  6.  Temperature  Profiles  at  £  ■  0,  0.1.  0.2,  and  0.3 


Figure  6  also  shows  that  because  of  heat  addition  to  the  fluid  from  the 
wall,  the  surface  temperature  for  the  Isothermal  wall  is  higher  than  that  for 
the  adiabatic  wall  by  3.34  percent  at  5  -  0.1,  6.6  percent  at  C  ■  0.2,  and  10.9 
percent  at  (  •  0.3;  giving  the  temperature  ratio  Tw  (isothermal)  to  Tw 
(adiabatic)  of  approximately  (1  +  0.3350*  Hall  considered  the  case  with  a 
finite  wall  thermal  conductivity.^  He  found  that  the  temperature  in  the 
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boundary  layer  on  an  adiabatic  wall  Is  the  same  as  the  temperature  in  the 
lnvlscid  freestreaa.  This  result  is  equivalent  to  the  zero-order  solution 
(l.e.,  £  ■  0)  of  Eq.  (40). 


C.  Heat  Transfer 


For  an  Isothermal  wall  the  local  heat  transfer  from  the  wall  to  the  fluid 
Is 


Qw 


.  3T 
3Y  Y-0 


(41) 


Written 


in  dimensionless  form. 


36 

3y'  y’-O 


it  becomes 


(42) 


Transformed  into  (£,z)  coordinate  system,  it  is 

qw  -  -t"1/2yH5~1/2(G  -  [g;(0)+5gj(0)+€2g£(0) ]  (43) 

The  values  of  g^,  g| ,  and  g£  are,  respectively,  0.9575,  0.15074,  and  -0.25445. 
Figure  7  and  Table  3  show  that  at  a  given  instant  the  heat  transfer  increases 
rapidly  behind  the  wavefront  until  it  approaches  a  maximum,  and  then  likely 
decreases.  This  can  be  explained  as  a  result  of  an  interplay  among  the 
expansion  of  the  freestream,  the  growing  boundary-layer  thickness,  and  the 
viscous  heating.  From  Eq.  (2)  we  notice  that  the  freestream  temperature 
decreases  with  increasing  distance  from  the  wavefront.  Consequently,  heat 
transfer  Is  from  the  wall  to  the  fluid.  On  the  other  hand,  the  combined  effect 
of  the  growing  boundary-layer  thickness  and  the  viscous  heating  provides  a 
resistance  to  the  heat  transfer  further  downstream.  Equation  (43)  shows  that 
the  heat  transfer  at  a  given  location  decreases  with  respect  to  time  and  is 
approximately  proportionate  to  t-*^.  Physically  this  is  reasonable  because  at 
a  given  location  the  boundary  layer  grows  in  time  after  the  expansion  wave 
passes  by.  The  growth  of  the  boundary  layer  deters  the  heat  transfer. 


Figure  7.  Skin  Friction  and  Heat  Transfer 


TABLE  3.  SKIN  FRICTION  AND  HEAT  TRANSFER 


C.  (Jtcf) 

o.r" 

1 

bothcnnal 

Adiabatic 

Isothermal 

0 

0.0000 

0.0000 

0.0000 

0.023 

0.1409 

0.1403 

0.0491 

0.030 

0.2034 

0.2042 

0.0073 

0.073 

0.2407 

0.2403 

0.0003 

0.100 

0.2039 

0.2003 

0.0901 

0.123 

0.3130 

0.3091 

o.om 

0.130 

0.3390 

0.3337 

0.1041 

0.173 

0.3027 

0.3353 

0.1090 

0.200 

0.3033 

0.3743 

0.1130 

0.223 

0.4010 

0.3912 

0.1101 

0.230 

0.4103 

0.4003 

0.1100 

0.273 

0.4337 

0.4199 

0.1203 

0.300 

0.4470 

0.4321 

0.1210 

0.323 

0.4003 

0.4431 

0.1220 

0.375 

0.4024 

0.4010 

0.1231 

0.400 

0.4920 

0.4097 

0.1220 

0.430 

0.3007 

0.4031 

0.1214 

0.300 

0.5224 

0.4939 

0.1191 

The  coefficient  of  local  heat  transfer  is  determined  by  Newton's  cooling 

law 


0,  -  MW 
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h  -  qwR(Re  k)~A/Z[ (1/y)-g]~ 


(44) 


where  R  Is  the  gas  constant.  The  function  G  Is  given  In  Eq.  (2). 
D.  Skin  Friction 

The  friction  coefficient  on  the  wall  is  defined  as 


_1 _  -  3U 

Cf  “  n  2  W  3Y  j  Y-0 
pU 


(45) 


After  coordinate  transformation  and  series  expansion,  Eq.  (45)  becomes 


Cf  -  (Re  tO”1/2  FH[  fg(0)  +  Cfj'(O)  +  ?2f“(0)]  (46) 

The  values  of  fQ,  fj'(O),  and  f£(0)  are,  respectively,  1.1284,  0.7946,  0.41781 
for  the  isothermal  wall  and  1.1284,  0.65896,  0.32903  for  the  adiabatic  wall. 

The  plot  of  Cj(Ret)1^  is  also  given  in  Figure  7  and  Table  3.  There  is  no  large 
difference  in  skin  friction  between  the  isothermal  wall  and  the  adiabatic  wall 
because  both  have  very  similar  velocity  profiles  as  shown  in  Figures  4  and  5. 
Specifically,  at  C  ■  0.3  the  skin  friction  for  the  isothermal  wall  is 
approximately  5.5  percent  higher  than  that  for  the  adiabatic  wall. 


V.  SUMMARY  AND  CONCLUSIONS 

Solutions  were  obtained  for  the  unsteady  compressible  laminar  boundary- 
layer  flow  that  develops  within  a  centered  expansion  wave  for  isothermal  and 
adiabatic  walls.  The  solutions  were  obtained  by  a  method  of  similarity 
transformation  via  one-parameter  groups  and  a  power  series  expansion  in  terms 
of  a  dimensionless  distance  £  from  the  expansion  wavefront.  The  series 
expansion  includes  terms  up  to  second  order  of  €,  and  thus  constitutes  an 
improvement  over  Hall's  zero-order  solution  for  the  temperature  in  the  boundary 
layer  over  an  adiabatic  wall.  The  present  results  show  that  at  a  given 
location  the  isothermal  wall  surface  temperature  is  higher  than  the  adiabatic 
wall  surface  temperature  by  a  factor  of  approximately  0.335  £• 
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Function  F*  of  Eq.  (33) 


On 


NOMENCLATURE 


coefficient  of  local  skin  friction 
specific  heat  at  constant  pressure 
invi8cid  velocity  function,  Eq.  (1) 
right-hand  side  of  Eq.(33),  Appendix  A 
transformed  streamf unction,  Eq.  (21) 
transformed  streamfunctlon,  Eq.  (26),  n  *  1, 
invlscld  temperature  function,  Eq.  (2) 
right-hand  side  of  Eq.  (34),  Appendix  A 
transformed  temperature  function,  Eq.  (22) 
transformed  temperature  function,  Eq.  (28) 
invlscld  pressure  function,  Eq.  (3) 
heat  transfer  coefficient 
thermal  conductivity 
reference  length 
integer,  n  -  0,1,2,... 
dimensionless  pressure 
dimensional  pressure 
dimensionless  freestream  pressure 
Prandtl  number 
Reynolds  number 

dimensional  heat  flux  at  the  wall 
dimensionless  heat  flux  at  the  wall 


gas  constant 


dimensional  temperature 


“  dimensional  freestream  temperature 

■  dimensional  temperature  of  the  undisturbed  fluid 

*  reference  temperature  Ur  /R  ■  y  TQ 

-  dimensional  wall  temperature 

-  dimensionless  time  tll/L 

r 

■  dimensional  time 

-  dimensionless  time.  Figure  1 

■  x-component  velocity,  dimensional 

*  reference  velocity,  speed  of  the  expansion  wavefront  /yRTQ 

*  x-component  velocity,  dimensionless,  Eq.  (6) 

■  dimensionless  freestream  velocity 

-  y-component  velocity,  dimensional 

-  y-component  velocity,  dimensionless,  Eq.  (6) 

■  dimensionless  coordinates.  Figure  1  and  Eq.  (6) 

*  dimensional  coordinates 

■  dimensionless,  transformed  coordinate  y',  Eq.  (11) 

*  variable,  Eq.  (25) 

■  ratio  of  specific  heats 

■  similarity  variable,  Eq.  (20) 

*  dimensionless  temperature 

■  dimensionless  freestream  temperature 

*  dimensionless  wall  temperature 


dimensionless  viscosity 


•O  I 


w 


dimensional  viscosity 


■  reference  viscosity 
£  “  similarity  variable,  Eq.  (4) 

p  -  dimensionless  density,  Eq.  (6) 

-  dimensional  density 
p  ^  -  reference  density 

^  *  streamf unction,  Eq.  (12) 

•  *  outer  edge  of  the  boundary  layer 

Subscripts 
1  *  lnvlscld 

r  •  reference 

w  ■  wall  condition 


0 


»  undisturbed 
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